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Abstract 

Computer  generated  2D  microstructures  of  varying  second-phase  area  fraction 
(5%-30%),  aspect  ratio  (1-16),  and  degree  of  alignment  (where  the  reinforcement 
major-axis  orientation  is  random,  perfectly  aligned,  or  semi-aligned)  are  analyzed 
via  the  isotropic  and  directional  forms  of  the  computationally  efficient  Multi-Scale 
Analysis  of  Area  Fractions  (MSAAF)  technique.  The  impact  of  these  microstructure 
parameters  on  the  representative  volume  element  (RVE)  necessary  to  characterize 
a  microstructure  is  ascertained  with  variations  in  isotropic  and  directional  homoge¬ 
nous  length  scales,  derivative  quantities  of  the  MSAAF  technique.  Analysis  of  these 
results  produces  empirical  expressions  for  the  directional  homogenous  length  scale 
as  a  function  of  area  fraction  and  aspect  ratio  for  the  limiting  cases  of  random 
and  “perfect”  second  phase  alignment.  Generally,  particle  alignment  is  observed  to 
increase  the  aspect  ratio  of  a  microstructure’s  RVE  -  a  trend  amplified  by  higher  re¬ 
inforcement  aspect  ratios  and  lower  area  fractions.  Particle  alignment  also  decreases 
the  absolute  size  of  such  an  element  by  reducing  the  directional  homogenous  length 
scales  transverse  to  the  axis  of  alignment.  Periodic  boundary  conditions  on  the 
perimeter  of  the  synthetic  microstructures  are  used  to  characterize  the  error  in  the 
MSAAF  technique  via  multiple  instantiations  of  the  same  microstructure,  which 
further  indicates  that  the  statistical  variation  in  linear  fraction  measured  using  the 
directional  technique  error  in  the  directional  homogenous  length  scale  can  be  an 
order  of  magnitude  less  than  that  measured  by  the  isotropic  technique  homogenous 
length  scale  at  the  same  length  scale. 
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1  Introduction 


When  assessing  structure-property  relationships  for  various  materials,  it  is 
necessary  to  define  a  representative  length  scale  or  volume  element  for  charac¬ 
terization  or  simulation  to  avoid  misleading  predictions  of  macroscopic  defor¬ 
mation,  fracture,  or  transport  behavior.  Likewise,  a  model  or  material  sample 
must  have  sufficient  resolution  to  capture  the  phenomena  of  interest.  Ideally,  a 
component-sized  material  model  or  sample  should  have  the  highest  resolution 
possible;  however,  practical  limitations  on  experimental  and  computational 
resources  bound  the  extent  and  resolution  of  material  models. 

The  balance  between  these  two  often  competing  needs  of  extent  and  resolution 
has  been  an  underlying  theme  of  numerous  studies  regarding  characterization 
of  heterogeneous  materials.  Researchers  have  focused  on  developing  techniques 
for  determining  representative  length  scales  for  the  extent  of  microstructures 
based  on  area  fraction  or  volume  fraction  (1;  2;  3;  4);  such  approaches  have 
shown  sensitivity  to  microstructure  clustering  and,  in  fact,  are  excellent  met¬ 
rics  for  characterizing  this  property.  The  inverse  problem  of  segmenting  a  mi¬ 
crostructure  to  varying  length  scales  (homogenizing)  has  also  been  pursued  as 
a  route  for  identifying  the  critical  length  scale  at  which  macroscopic  properties 
are  controlled  by  shorter  length  scale  features;  these  features  are  often  associ¬ 
ated  with  second  phase  clustering,  its  absence,  or  contribution  to  percolation 
(5;  6;  7;  8).  In  fact,  recent  advances  in  multiscale  modeling  have  emerged  due 
to  the  ability  to  homogenize  some  regions  while  considering  the  local  behavior 
in  others;  such  techniques  incorporate  various  length  scales  associated  with 
the  microstructure  and  its  local  properties  in  their  formulation  (9;  10). 

The  correlation  between  critical  length  scales  in  the  microstructure  and  the 
inherent  mechanical  behavior  of  heterogeneous  microstructures  has  been  ex¬ 
amined  in  numerous  computational  and  experimental  studies.  For  example, 
Spowart  (11)  predicted  that  clustering  of  reinforcement  particles  in  an  Al- 
27.5%SiC  metal  matrix  composite  has  a  noticeable  effect  on  the  yield  strength 
and  an  even  stronger  effect  on  strain-hardening,  even  though  the  same  level  of 
clustering  has  only  a  weak  effect  on  the  elastic  modulus.  These  results  indicate 
that  each  aspect  of  mechanical  behavior  is  controlled  by  a  different  representa¬ 
tive  length  scale  in  the  microstructure  and,  consequently,  that  representative 
volume  elements  (RVEs)  of  increasing  size  were  necessary  to  characterize  elas¬ 
tic  modulus,  yield  strength,  and  strain-hardening,  respectively.  Additionally, 
Borbely  et  al.  (2)  measured  the  local  volume  fractions  in  an  Al-20%Al2O3 
metal-matrix  composite  through  3D  microtomography  to  define  a  microstruc¬ 
ture  correlation  length  and,  hence,  a  geometric  RVE  for  this  material.  Further 
simulations  of  elastic  and  plastic  behavior  with  various  window  sizes  showed 
that  the  RVE  necessary  to  obtain  accurate  effective  plastic  properties  was 
approximately  twice  the  size  of  the  RVE  needed  to  obtain  equally  accurate  ef- 
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fective  elastic  properties.  Moreover,  the  influence  of  second  phase  distribution 
on  damage  accumulation  and  fracture  has  also  been  characterized  by  numer¬ 
ous  authors  in  various  materials,  e.g.,  (12;  13;  14).  In  their  seminal  work  on 
crack  deflection  around  second  phase  inclusions  in  a  brittle  matrix,  Faber  and 
Evans  analyzed  and  experimentally  verified  that  a  clustered  second  phase  dis¬ 
tribution  -  as  defined  by  a  divergence  from  a  uniform  distribution  for  a  given 
volume  fraction  -  significantly  increased  fracture  resistance  in  ceramic  matrix 
composites  by  increasing  the  degree  of  crack  deflection  and  crack  twist  (12;  13). 
The  clustering  of  particles  is  associated  with  a  critical  length  scale  for  fracture. 
In  experiments  relating  to  ductile  fracture,  Wilks  analyzed  second  phase  in¬ 
clusions  in  a  ductile  matrix  and  found  that  deformation  processing  decreased 
the  size  of  the  RVE  in  an  Al-25%SiC  composite;  the  smaller  RVE  correlated 
with  (i)  an  increase  in  interparticle  separation  due  to  cluster  breakdown,  (ii) 
a  larger  length  scale  with  respect  to  the  ductile  fracture  process,  and  (iii)  a 
substantial  increase  in  fracture  toughness  (14).  Finite  element  simulations  of 
microstructure  domains  reconstructed  from  experimental  images  of  multiphase 
materials  have  also  provided  insight  into  plasticity,  damage  evolution,  and  the 
fracture  process  (15;  16;  17;  18;  19).  In  all  of  these  experiments  and  simula¬ 
tions,  the  length  scale  of  the  underlying  microstructure  plays  a  commanding 
role  in  the  mechanical  behavior  of  the  material. 


Therefore,  it  is  important  to  have  a  consistent  methodology  for  characterizing 
the  critical  length  scales  within  heterogeneous  materials.  The  focus  of  this  work 
is  the  effect  of  second-phase  anisotropy  on  critical  length  scales  and  the  cor¬ 
responding  RVE  size.  Orthotropic  (aligned)  microstructures  are  very  common 
in  materials  where  deformation  processing  has  been  imposed  or  directional 
synthesis  routes  have  been  pursued  (14).  These  microstructures  constitute  a 
sound  basis  for  analyzing  limiting  bounds  on  anisotropic  length  scale  effects. 
Moreover,  such  microstructures  represent  a  basis  for  analyzing  microstructures 
with  several  preferred  directions  of  alignment  (e.g.,  precipitated  microstruc¬ 
tures  with  preferred  orientation  relationships) ,  or  more  generalized  anisotropy. 
In  this  work,  we  investigate  the  interaction  between  length-scale  effects  and 
second  phase  orthotropy  through  applying  the  isotropic  and  directional  forms 
of  the  Multi-Scale  Analysis  of  Area  Fractions  (MSAAF)  technique  (1;  20) 
to  computer-generated  two-dimensional  microstructure  images  containing  el¬ 
lipsoidal  particles  of  varying  aspect  ratio,  area  fraction,  and  propensity  for 
alignment. 
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Table  1 

The  test  matrix  for  these  simulations. 


Factor 

Factor  Levels 

Volume  Fraction 

5% 

10% 

20% 

30% 

Aspect  Ratio 

1 

2 

4 

8 

16 

Orientation  Distribution 

90° 

90°  ±  20° 

Random 

2  Methodology 

2.1  Synthetic  micro  structure  conditions 


Synthetic  microstructures  were  generated  with  the  factors  and  factor  levels 
summarized  in  Table  1.  The  four  area  fractions  (Af)  and  five  aspect  ratios 
(AR)  were  chosen  to  give  a  wide  range  and  sufficiently  resolved  gradient  in 
microstructure  characteristics  for  the  subsequent  analyses.  The  orientation  dis¬ 
tributions  studied  include  limiting  conditions  where  all  particles  are  aligned 
vertically  (90°),  all  particles  are  randomly  oriented  (Random),  and  an  inter¬ 
mediary  factor  level  (90°  ±20°)  that  reflects  orientations  normally  distributed 
about  90°  with  a  standard  deviation  of  20°.  In  all,  60  distinct  synthetic  mi¬ 
crostructures  (conditions)  were  studied  corresponding  to  the  permutation  of 
these  factors. 


2.2  Synthetic  microstructure  generation 


Each  synthetic  microstructure  condition  was  generated  via  a  simple  proce¬ 
dure,  similar  to  the  Random  Sequential  Adsorption  (RSA)  technique,  e.g., 
(21;  22).  First,  particle  dimensions  (in  pixels 1 )  were  calculated  for  each  aspect 
ratio  along  with  the  number  of  particles  necessary  to  obtain  the  prescribed 
area  fraction.  In  order  to  mitigate  any  unintended  effect  of  particle  size  dis¬ 
tribution,  particle  size  (area)  was  kept  constant  through  successively  higher 
area  fractions;  the  details  of  this  effect  will  be  considered  in  future  work  as  a 
subset  of  the  more  general  interaction  of  microstructure  anisotropy  and  par¬ 
ticle  polydispersivity.  Particle  orientation  was  then  randomly  sampled  from 
the  predetermined  orientation  distribution  while  the  coordinates  of  the  par¬ 
ticle  centroid  were  selected  via  a  random  number  generator.  The  influence  of 
boundary  effects  on  particle  placement  was  mitigated  by  imposing  periodic 

1  In  these  synthetic  microstructures,  absolute  length  scale  is  arbitrary;  since  fea¬ 
tures  dimensioned  in  pixels  may  seem  unphysical,  the  reader  may  wish  to  consider 
a  scale  of  1  pixel  =  1  yrn. 
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boundary  conditions  on  the  vertical  and  horizontal  edges  of  the  microstruc¬ 
ture  domain;  particles  that  extend  beyond  the  image  bounds  were  wrapped 
back  onto  the  opposite  side  of  the  image.  Analytic  criteria  were  used  to  re¬ 
ject  overlapping  particles.  Particle  placement  proceeded  until  the  specified 
area  fraction  was  obtained.  Resulting  microstructures  contained  a  minimum 
of  4,000  particles  and  were  of  4,096  x  4,096  pixels  in  extent. 

Sample  regions  from  each  synthetic  microstructure  condition  as  a  function  of 
AR  and  Af  are  depicted  in  Fig.  1.  The  general  influence  of  each  considered 
orientation  distribution  on  a  high  area  fraction  (Af  =  30%),  high  aspect  ratio 
(AR  =  16)  microstructure  is  shown  in  Fig.  2.  Notice  in  these  two  figures  that 
particles  with  similar  orientations  cluster  at  higher  area  fractions  and  aspect 
ratios.  This  effect  can  be  correlated  with  the  variation  in  jamming  threshold 
(23)  for  such  microstructures  and  indicates  that  the  distribution  of  particle 
centroids  may  not  be  truly  random,  due  to  the  stochastic  nature  of  particle 
placement. 


2.3  Multi-Scale  Characterization 

2.3.1  Isotropic  Microstructures 

The  homogenous  length  scale  (LH)  was  developed  to  determine  a  representa¬ 
tive  length  scale  for  a  particular  microstructure  based  on  the  variations  in  area 
fraction  of  a  second  phase  as  a  function  of  length  scale  (Q)  (1).  Qualitatively, 
Lh  for  a  particular  microstructure  is  the  length  scale  at  which  the  variation 
between  each  sub-region  becomes  statistically  indistinguishable  from  a  larger 
area  of  material.  In  other  words,  Lh  is  the  minimum  length  scale  necessary 
to  construct  a  volume  element  representative  of  a  particular  microstructure  to 
a  given  degree  of  statistical  confidence,  typically  selected  to  be  99.0%  in  the 
isotropic  case. 

The  homogeneous  length  scale  of  a  real  or  synthetic  microstructure  can  be 
quantified  by  a  technique  termed  the  multi-scale  analysis  of  area  fractions 
(MSAAF)  (1).  In  this  method,  a  digitized  microstructure  is  re-sampled  at 
various  resolutions  (corresponding  to  various  length  scales,  Qn)  to  determine 
the  area  fraction  of  a  particular  sub-region.  Take  for  example,  the  first  column 
of  images  in  Fig.  3,  where  an  initial  microstructure  is  subdivided  into  3  lower 
resolutions/length  scales:  Q i,  Q2,  Qz-  After  subdivision,  the  effective  area 
fraction  (Af)  is  assigned  to  each  sub-region,  visually  depicted  as  a  gray  level 
in  Fig.  3.  By  using  all  sub-regions  available  at  a  chosen  length  scale,  an  average 
area  fraction  and  standard  deviation  (<74  )  in  area  fraction  is  determined  for 
that  Qn.  Though  the  average  area  fraction  over  all  sub-regions  for  a  particular 
length  scale  will  not  change  with  length  scale,  the  standard  deviation  over  all 
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Fig.  1.  All  synthetic  microstructure  conditions  for  randomly  oriented  particles  displayed  as  a  function  of  aspect  ratio  ( AR )  and  area 
fraction  ( Af ). 


(a)  90°  (b)  90°  ±  20°  (c)  Random 

Fig.  2.  The  influence  of  orientation  distribution  on  synthetic  microstructures  con¬ 
taining  30%  area  fraction  second  phase  with  a  16:1  aspect  ratio.  The  512  x  512 
subimages  show  the  alignment  of  particles  in  the  (a)  90°,  (b)  90°  ±  20°,  and  (c) 
random  orientation  distributions. 

sub-regions  will. 

Consequently,  the  principal  result  of  this  technique  is  a  plot  depicting  the 
evolution  of  the  coefficient  of  variation  (fl/)  -  the  ratio  of  standard  deviation  in 
second-phase  area  fraction  to  the  global  area  fraction,  -  as  a  function 

of  Q.  Such  a  plot  for  the  20%A/-16AR  condition  is  depicted  in  Fig.  4.  When 
this  plot  decreases  to  a  specified  level  of  confidence  ( T ) ,  the  corresponding 
Q  is  termed  the  homogenous  length  scale  (Lm).  Differences  in  homogeneity 
between  isotropic  microstructures  can  be  characterized  using  a  value  of  LH 
for  each  microstructure  at  a  predetermined  level  of  confidence  (1).  For  this 
reason,  we  use  the  isotropic  MSAAF  analysis  of  synthetic  microstructures  to 
characterize  homogeneity  with  the  1%  homogenous  length  scale,  L °If1,  i.e., 
the  scale  at  which  the  variation  in  microstructure  sub-region  area  fraction 
becomes  less  than  1%  (4/  <  0.01). 


2.3.2  Anisotropic  Microstructures 

When  studying  a  area  of  material  that  is  both  anisotropic  and  inhomogeneous, 
it  is  useful  to  logically  extend  the  concept  of  the  homogenous  length  scale  by 
measuring  the  directional  variation  in  linear  area  fraction  as  a  function  of 
length  scale  to  compute  a  directional  homogenous  length  scale,  Lhh ,  where 
n  is  an  integer  index  (valued  from  1-3  for  a  3-dimensional  microstructure) 
denoting  a  principal  particular  (orthogonal)  direction  in  the  material  (20). 
Only  a  simple  modification  to  the  MSAAF  technique  is  required  to  obtain 
values  of  LHn  from  plane  sections  of  a  material.  As  shown  in  the  second  and 
third  columns  of  Fig.  3,  rather  than  subdividing  the  microstructure  into  square 
elements,  each  line  of  pixels  is  subdivided  into  a  number  of  strips  for  each 
corresponding  length  scale  of  interest  ( Qn ).  The  variation  in  area  fraction  of 
each  of  these  strips  is  used  to  generate  directional  MSAAF  plots  (also  depicted 
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(a)  Isotropic 


(b)  Vertical 


(c)  Horizontal 


Fig.  3.  Evolution  of  the  density  area  fraction  of  second  phase  (depicted  by  gray  level) 
as  a  function  of  length  scale  Q  according  to  the  (a)  isotropic  MSAAF  technique, 
as  well  as  the  directional  technique  applied  in  the  (b)  vertical  and  (c)  horizontal 
directions  for  the  synthesized  20% A f -16 AR  microstructure. 


Quilt  Dimension,  Q 
(a) 


(b) 

Fig.  4.  Isotropic  and  directional  MSAAF  plots  for  the  synthesized  20%Af-16AR 
microstructure  with  (a)  aligned  particles  and  (b)  randomly  oriented  particles. 


for  the  20%A/-16AR  in  Fig.  4)  to  determine  the  orthogonal  10%  (\F  =  0.1) 
homogenous  length  scales,  Lhi  (vertical)  and  Lh2  (horizontal).  A  value  of 
=  0.1  (he.,  90%  confidence)  is  chosen  due  to  practical  limitations  on  the 
sizes  of  images  that  can  typically  be  obtained  either  experimentally  or  by 
simulation. 
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For  all  synthetic  microstructure  conditions  resulting  from  the  permutations  of 
the  factors  in  Table  1,  the  MSAAF  technique  and  its  directional  counterpart 
were  used  to  assess  the  influence  of  these  factors  on  LH ,  LH 1}  and  LH2.  Since 
the  synthetic  microstructure  images  are  periodic,  multiple  instantiations  of  the 
same  microstructure  were  sampled  (by  re-centering  the  image)  to  characterize 
the  inherent  variation  in  the  isotropic  and  directional  MSAAF  techniques.  To 
this  end,  30  distinct  instantiations  of  the  same  synthetic  microstructure  were 
used  to  calculate  error  values  (3  standard  deviations)  at  each  length  scale 
considered  in  the  MSAAF  analysis,  as  is  also  depicted  by  the  error  bars  in 
Fig.  4. 


3  Results  and  Discussion 


3.1  Isotropic  and  Directional  MSAAF  Results 


The  isotropic  homogenous  length  scale,  Lh( 0.01),  for  different  area  fractions 
and  aspect  ratios  of  aligned  particles  (90°  orientation  distribution)  is  depicted 
in  Fig.  5.  First,  as  expected  from  previous  work  (14),  there  is  a  general  decrease 
in  Lh  with  increasing  second  phase  area  fraction.  Second,  there  is  seemingly 
no  monotonic  correlation  between  L#(0.01)  and  aspect  ratio.  Strongly  con¬ 
tributing  to  this  absence  of  correlation  is  the  broad  variation  in  Lh (0.01)  be¬ 
tween  different  instantiations  of  each  microstructure  that  often  overlaps  with 
adjacent  conditions. 

Directional  homogenous  length  scales  LH\( 0.1)  and  LH2( 0.1)  for  the  aligned 
microstructures  (90°  orientation  distributions)  are  depicted  in  Figs.  6(a)  and 
6(b),  respectively.  Here,  trends  with  both  area  fraction  and  aspect  ratio  are 
apparent.  Again,  the  inverse  relationship  between  homogenous  length  scale 
and  area  fraction  is  clear.  However,  the  directional  length  scale  along  the  axis 
of  alignment  (LH i)  increases  with  aspect  ratio  while  the  opposite  trend  is 
observed  in  the  transverse  ( Lh2 )  direction. 

Comparing  Figs.  5  and  6  clearly  indicates  that  for  the  isotropic  MSAAF  tech¬ 
nique,  error  in  Lh  measurement  can  be  significant  when  comparing  orthotropic 
microstructures.  This  is  due  in  part  to  the  (m  —  1)  character  of  the  standard 
deviation  in  area  fraction  used  to  compute  T.  where  m  is  the  number  of  sub- 
regions;  at  large  length  scales,  this  deviation  can  influence  the  measurement 
of  LH-  In  contrast,  the  variation  at  each  length  scale  resulting  from  the  direc¬ 
tional  MSAAF  technique  is  quite  small  since  even  at  the  largest  length  scale 
there  are  still  4096  sub-regions;  the  different  length  scales  more  accurately 
capture  the  contribution  of  orthotropy  in  the  microstructure.  For  this  reason, 
the  terms  containing  fewer  than  16  only  1  or  4  sub-regions  are  dropped  when 
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Particle  Aspect  Ratio,  AR 


Fig.  5.  Isotropic  homogenous  length  scale,  Lh(0.01),  as  a  function  of  second-phase 
aspect  ratio  and  area  fraction  for  aligned  (90°  orientation  distribution)  particles; 
error  bars  represent  the  3cr  deviation  from  30  instantiations  of  each  condition. 

calculating  the  isotropic  LH  to  avoid  including  spurious  values  of  LH. 

The  effect  of  orientation  distribution  on  the  directional  homogenous  length 
scales  as  a  function  of  aspect  ratio  for  the  synthetic  microstructures  containing 
30% A f  second  phase  is  captured  in  Fig.  7.  As  expected,  there  is  no  difference 
in  directional  length  scales  for  the  random  condition.  However,  the  increas¬ 
ing  propensity  for  second-phase  alignment  increases  the  separation  of  these 
orthogonal  length  scales  by  increasing  LHi  and  decreasing  LH2.  This  trend  is 
characteristic  of  the  other  second  phase  area  fractions  as  well. 


3.2  Model  for  Homogenous  Length  Scales  in  Orthotropic  Micro  structures 


MSAAF  data  has  previously  been  fit  (1)  to  an  equation  of  the  form 

1 

\M//(i  ~  A/0  +  a(2* 

where  a  is  a  geometric  factor  related  to  particle  size  and  area  fraction,  £ 
is  the  slope  of  the  long  length-scale  portion  of  the  MSAAF  plot,  and  Q*  is 
the  discretized  length  scale  (in  pixels)  (1).  In  order  to  understand  the  effect 
of  microstructure  factors  on  homogenous  length  scales,  least-squares  fits  to 
Eq.  1  were  used  to  characterize  the  behavior  of  a  and  £  for  each  synthetic 


11 


Particle  Aspect  Ratio,  AR 
(a) 


Particle  Aspect  Ratio,  AR 

(b) 

Fig.  6.  Directional  homogeneous  length  scales,  (a)  Lh i(0.1)  and  (b)  Lh 2(0.1),  as  a 
function  of  second-phase  aspect  ratio  and  area  fraction  for  aligned  (90°  orientation 
distribution)  particles. 
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Particle  Aspect  Ratio,  AR 

Fig.  7.  Directional  homogenous  length  scales,  Lhi(0.1)  and  L# 2(0.1),  as  a  function 
of  second-phase  aspect  ratio  for  random,  semi-aligned,  and  aligned  (90°  orientation 
distribution)  particles  in  a  30% 71/  microstructure. 

microstructure  condition.  During  such  fits,  it  was  observed  that  while  Q*  =  Q 
fits  MSAAF  curves  where  a  is  small,  large  a  values  require  Q*  =  Q  —  1.  This  is 
simply  due  to  the  discretized  nature  of  the  image,  i.e.,  the  lower  limit  on  Q  is  1 
pixel.  This  modification  reflects  that  at  a  length  scale  of  Q  =  1,  the  coefficient 
of  variation  is  entirely  dependent  on  second  phase  area  fraction  (i.e.,  T  ( Af ) 
=  \J  A  —  Af )  /Aj ) .  Use  of  this  transformation  for  any  a  fully  reconciles  the 
short  and  long  length-scale  MSAAF  behaviors  originally  proposed  by  Spowart 
et  al.  (1)  and  greatly  improves  the  fit  of  Eq.  1  to  MSAAF  results,  suggesting 
that  the  relationship 


V{Af)  =  7  (2) 

\! Af/ (1  -  A/)  +  a(Q  -  1)  2€ 

is  more  universal  for  discrete  data  than  Eq.  1,  particularly  at  short  length 
scales. 

The  solid  lines  for  the  directional  MSAAF  plots  shown  in  Fig.  4  are  typical 
of  fitting  Eq.  2  to  MSAAF  results  for  all  synthetic  microstructures.  Over  all 
conditions  isotropic  ~  —  1  and  £ directional  ~  —0.5,  indicating  that  £  is  insensitive 
to  all  considered  microstructure  factors,  as  well  as  the  method  (isotropic  or 
directional)  used  for  characterizing  homogeneity.  On  the  other  hand,  a  var¬ 
ied  significantly  with  second  phase  area  fraction,  aspect  ratio,  and  degree  of 
alignment.  Moreover,  it  has  previously  been  shown  (1)  (11)  that  £  can  vary 
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Table  2 

Coefficients  for  the  a  model  of  Eq.  3  for  each  Directional  homogenous  length  scales, 
(a)  Lhi,  and  (b)  Lh2  as  a  function  of  second-phase  aspect  ratio  and  area  fraction 
for  aligned  (90°  orientation  distribution)  particles. 


Orientation 

co 

Cl 

C2 

C3 

C4 

Random,  Lhi 

0.0201 

0.2491 

0.0117 

0.0056 

0.0327 

Random,  Lh2 

0.0190 

0.2682 

0.0096 

0.0058 

0.0303 

90°  d=  20°  Aligned,  Lhi 

-0.0049 

0.1613 

-0.0010 

0.0024 

0.0733 

90°  d=  20°  Transverse,  Lh2 

0.0512 

0.2267 

0.0569 

0.0090 

-0.0066 

90°  Aligned,  Lhi 

-0.0091 

0.0983 

-0.0046 

0.0006 

0.0895 

90°  Transverse,  Lh2 

0.0583 

0.2488 

0.0570 

0.0119 

-0.0183 

significantly  with  degree  of  clustering,  although  this  was  not  studied  in  the 
present  effort.  Fig.  8  shows  the  evolution  of  a  as  function  of  second  phase 
area  fraction  and  aspect  ratio  for  the  (a)  random  orientation  distribution,  and 
the  “perfectly”  aligned  orientation  distribution  in  the  (b)  LHi  and  (c)  LH2 
directions. 

Inspecting  the  results  for  the  random,  partially  aligned  (90°  ±20°),  and  “per¬ 
fectly”  aligned  (90°)  conditions  showed  that  a  obeys  a  relationship  of  the 
form 


a  —  Af 


Co  +  c\Af  +  c2AfAR  +  c$AR  +  C4 


(3) 


where  cn  are  constants  that  are  tabulated  in  Table  2  for  each  orientation  dis¬ 
tribution  considered  in  this  study.  The  quality  of  the  fit  from  these  coefficients 
to  the  calculated  a  values  is  depicted  in  Fig.  8(a)-(c)  by  the  solid  lines. 


3.3  Implications  for  Representative  Volume  Elements  Size  in  Orthotropic 
Materials 


If  it  is  assumed  that  £  =  —0.5,  then  Eq.  2  can  be  rearranged  such  that  LH, 
Lm,  and  Lh2  can  be  characterized  by  a  relationship  of  the  form 


1 


Af 

1  -  Af_ 


a 


1 


Af 

1  -  Af_ 


(4) 


14 


a  for  Random  Alignment 


Particle  Aspect  Ratio,  AR 

Particle  Aspect  Ratio,  AR 

Particle  Aspect  Ratio,  AR 

(a) 

(b) 

(c) 

Fig.  8.  Evolution  of  the  parameter  a  as  a  function  of  second  phase  area  fraction  and  aspect  ratio  for  (a)  random  and  “perfectly”  aligned 
(90°)  orientation  distributions  in  the  (b)  Lin(O.l)  and  (c)  ^if2(0.1)  directions. 


where  the  index  n  of  Eq.  4  is  1,  2,  or  3  depending  on  the  principal  direction  in 
the  microstructure.  Such  an  expression  can  then  be  used  to  estimate  the  size  of 
a  representative  volume  element  for  the  considered  orthotropic  microstructures 
as  a  function  of  the  second  phase  volume  fraction  (Vf),  aspect  ratio,  and 
desired  level  of  confidence  (T),  i.e., 


v  RVE  ~  ~  „  „  „ 

QL1QL2QL3 


1 

\|/2 


1  -  V>J 


(5) 


If  symmetry  is  maintained  about  the  axis  of  alignment  (the  1  direction)  and 
the  second  phase  is  needle-shaped  not  disc- shaped,  then  a2  =  «3  and  the 
size  of  the  RVE  necessary  to  characterize  the  3-dimensional  variants  of  the 
orthotropic  microstructures  can  be  determined  as  a  function  of  the  factors 
considered  in  this  work.  To  that  end,  Fig.  9(a-c)  captures  the  effect  of  second 
phase  volume  fraction  and  aspect  ratio  on  RVE  size  for  the  (a)  random,  (b) 
90°  ±20°,  and  (c)  “perfectly”  (90°)  aligned  orientation  distributions. 


4  Discussion 


There  are  a  few  caveats  of  the  2D  synthetic  microstructures  generated  in 
this  work.  The  random  placement  of  particles  may  not  be  entirely  repre¬ 
sentative  of  extrusion  or  directional  synthesis  processing.  For  example,  the 
nearest  neighbor  separation  could  contract  during  extrusion  in  addition  to 
second  phase  alignment  or  there  could  be  alignment  of  second  phase  cen¬ 
troids  in  addition  to  the  orientation  of  the  principal  axis.  In  this  respect,  the 
methodology  for  generating  synthetic  microstructures  could  be  altered  to  bet¬ 
ter  represent  the  microstructure  within  a  specific  material  processed  in  a  spe¬ 
cific  manner,  i.e.,  two-point  correlation  functions  of  experimental  microstruc¬ 
tures  could  be  used  to  generate  statistically-similar  synthetic  microstructures 
(24;  25;  26;  27;  28;  29;  30).  Additionally,  there  are  many  aspects  of  microstruc¬ 
tures  in  real  materials  that  have  not  been  explored  in  the  current  work,  i.e., 
particle  morphology  (in  addition  to  aspect  ratio),  particle  size  distribution, 
and  particle  clustering  (spatial  distribution),  to  name  a  few.  We  leave  these 
microstructure-specific  aspects  of  second  phase  particles  for  future  work. 

Nevertheless,  by  controlling  the  second-phase  particle  alignment,  aspect  ratio, 
and  area  fraction  when  generating  the  2D  periodic  synthetic  microstructures, 
the  following  work  illuminates  how  the  computationally-efficient  MSAAF  tech¬ 
nique  can  be  used  for  understanding  the  inherent  alignment  of  second  phase 
particles  in  experimental  images.  The  two  microstructure  parameters,  a  and 
£,  used  for  fitting  Eq.  2  represent  the  evolution  of  the  variation  in  area  fraction 
as  a  function  of  length  scale.  The  directional  and  homogeneous  length  scales 
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(a)  (b)  (c) 

Fig.  9.  Effect  of  second  phase  area  volume  fraction  (Af)  and  aspect  ratio  (AR)  on  2-dimensional  RVE  size  for  the  (a)  random 
90°  ±20°,  and  (c)  “perfectly”  (90°)  aligned  orientation  distributions. 


from  the  MSAAF  technique  signify  a  quantitative  metric  for  characterizing 
the  representative  length  scale  or  volume  element  for  a  given  microstructure. 
The  significance  of  this  work  is  that  the  analysis  of  experimental  images  with 
the  MSAAF  technique  can  be  compared  and  interpreted  based  on  the  results 
of  various  synthetic  microstructures  examined  in  this  work. 

Future  work  will  encompass  generating  3D  periodic  synthetic  microstructures 
with  second  phase  particles  of  known  sizes,  aspect  ratios,  and  orientations;  an¬ 
alyzing  2D  orthogonal  slices  may  help  our  understanding  of  how  quantitative 
information  informed  from  2D  orthogonal  slices  relates  to  the  known  3D  struc¬ 
ture.  This  will  test  the  validity  of  stereological  assumptions  often  employed 
for  relating  2D  microstructural  information  to  3D  volumes  (31).  Additional 
work  will  concentrate  on  generating  synthetic  microstructures  with  various 
degrees  of  particle  clustering  to  improve  our  understanding  of  how  cluster¬ 
ing  impacts  the  representative  length  scales  of  realistic  microstructures.  Also, 
as  shown  in  the  present  work,  the  directional  length  scale  indicates  align¬ 
ment  (or  non-alignment)  in  a  particular  direction;  by  applying  the  directional 
MSAAF  technique  along  multiple  orientations,  the  associated  alignment  tex¬ 
ture  of  second-phase  particles  within  a  microstructure  image  can  be  easily 
ascertained. 


5  Summary 


In  this  work,  synthetic  two-phase  microstructures  with  known  area  fractions, 
aspect  ratios,  and  orientations  of  second  phase  particles  were  generated  to 
probe  the  effects  of  orthotropy  on  metrics  used  to  measure  representative 
length  scales:  the  isotropic  and  directional  homogenous  length  scales  calcu¬ 
lated  from  the  Multi-Scale  Analysis  of  Area  Fractions  (MSAAF)  technique. 
The  prime  result  of  this  work  is  a  template  for  determining  the  representative 
length  scale  (and  volume)  of  microstructure  in  a  given  2D  section  as  a  func¬ 
tion  of  specific  microstructure  factors  (area  fraction,  aspect  ratio,  alignment 
propensity).  Ancillary  results  of  this  work  are: 

(1)  The  calculated  standard  deviation  (error)  from  30  instantiations  of  each 
periodic  synthetic  microstructure  is  substantially  lower  for  directional  ho¬ 
mogenous  length  scales  when  compared  with  the  isotropic  homogenous 
length  scale.  This  is  attributable  to  the  number  of  constituent  subdivi¬ 
sions  of  an  image  at  larger  length  scales.  Therefore,  to  avoid  this  deviation 
the  extent  of  the  domain  must  be  large  compared  to  the  representative 
length  scale  being  measured;  this  is  especially  important  for  the  isotropic 
MSAAF  technique. 

(2)  The  directional  MSAAF  technique  is  useful  for  characterizing  the  rela¬ 
tive  alignment  of  second  phase  particles  in  the  synthetic  microstructures 
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generated.  The  deviation  between  the  directional  length  scales  in  two  or¬ 
thogonal  directions  increases  as  a  function  of  particle  alignment  and  as¬ 
pect  ratio.  The  homogeneous  and  directional  length  scales  increase  with 
decreasing  area  fraction  of  particles. 

The  present  work  shows  how  the  isotropic  and  directional  MSAAF  techniques 
can  help  quantitatively  characterize  orthotropic  microstructures  with  second 
phase  particles  of  varying  alignments. 
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